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,__! I Abstract 

\^ • Beginning with a brief sketch of the derivation of Hawking's theorem of horizon 

CN . area increase, based on the Raychaudhuri equation, we go on to discuss the issue 

^-p I as to whether generic black holes, undergoing Hawking radiation, can ever remain 

f~^ ■ in stable thermal equilibrium with that radiation. We derive a universal criterion 

^O . for such a stability, which relates the black hole mass and microcanonical entropy, 

^^ I both of which are well-defined within the context of the Isolated Horizon, and in 

^^' principle calculable within Loop Quantum Gravity. The criterion is argued to hold 

I . even when thermal fluctuations of electric charge are considered, within a grand 

{^1 canonical ensemble. 

>: 

X ■ 1 Introduction 

... the equation for the rate of change of expansion (Raychaudhuri's equation) 
plays a central role in the proofs of the singularity theorems. [J 

Even in its fifty-first year, the Raychaudhuri equation [2] retains its prime position as a 
tool to analyze spacetime singularities (for a review for non-specialists, see j3]). However, 
lest one should feel that the equation has only one specialized role in general relativity, 
namely that of delineating singularities, we remark that it has had far wider applications, 
as for example in establishing the so-called laws of black hole mechanics. The earliest of 
these laws or theorems is called Hawking's theorem of increase of horizon area of a black 
hole. It states that |lj 

The area of the event horizon of an isolated stationary black hole can never 
decrease in any physical process 
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Figure 1: Black hole spacetime emerging from spherically symmetric collapse of a massive 
star M; ellipses denote successive sizes of the collapsing star, and the null surface H is the 
event horizon 



We recall first that a black hole spacetime B is the set of events that lie in the comple- 
ment of the chronological past of events infinitely distant and infinitely far in the future 
(the future asymptotic null infinity X"*"). The event horizon of the black hole is the bound- 
ary dB of events in spacetime accessible to such observers. It is a null 2+1 dimensional 
surface hiding the black hole singularity. 

A sketch of the proof of Hawking's theorem may be given as follows [4J : Recall that 
the Raychaudhuri equation j2] 



dr 
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relates the evolution of the (volume) expansion 9{t) of (timelike and null) geodesies, to 
the 'shear' a which quantifies relative stretching of geodesies in a congruence, the rotation 
or vorticity u which exhibits how the geodesies twist around each other as they evolve 
and the spacetime curvature responsible for geodesic deviation. 

Consider any null surface generated by null geodesies which we assume complete. A 
fixed time slice of such a surface is a spacelike 2-sphere of area A. The Raychaudhuri 
equation ((T]) now implies that 
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The expansion 6 = dlogA/dr. Assume ^(0) > in contrast to the case for studying singu- 
larities; the area A then increases locally. Consider fixed time slices 5*1, 5*2 of event horizon 
at Ti,T2 > Ti. The number of null geodesic generators intersecting 5*2 is at least the same or 



bigger than those intersecting Si. At any instant, A oc no. of Null geodesic generators; 
it follows then that =^ A2 > Ai which is what the theorem states. 

Hawking's theorem brings to mind the law of entropy increase (for isolated systems) 
that follows from the second law of thermodynamics. Yet, for a classical black hole 
spacetime, the complete absence of 'microstates' involving atoms, electrons, photons or 
other discrete entities, stymies pushing the analogy very far. The situation turns curiouser 
because of other 'laws of (stationary) black hole mechanics' derived on the basis of general 
relativity |1], wherein the geometrical quantity known as surface gravity (k) appears as 
an analogue of temperature. 

This thermodynamic analogy of the theorems on stationary black hole mechanics, 
especially, Hawking's theorem, led Bekenstein [7j to postulate that black holes must have 
an entropy proportional to the area of their horizons 
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Sbh — C ^B -J2 ■ ('^) 



Here, C is a dimensionless constant of 0(1), and Ip is the Planck length Ip = [Gh/c'Y^'^ ~ 
10"'^'^ cm, characterizing the length scale at which spacetime can no longer be thought 
of in terms of classical Riemannian geometry. It follows that black hole thermodynamics 
must have to do with quantum, rather than classical, general relativity. The microstates 
from which S}yh originates must then correspond to the states of the quantum spacetime 
associated with the black hole, especially the event horizon. 

Similarly, black hole temperature also has quantum origins, being associated with 
the mysterious Hawking radiation [H]. Every generic black hole radiates particles in a 
thermal (black body) spectrum, at a temperature T{,^ = huhorizon, a surprising result, 
if one remembers that the surface gravity k is purely a geometric quantity. Hawking 
radiation thus corresponds to transformation of quantum spacetime states describing the 
black hole, into thermal particle states. In a sense, it is an inversion of gravitational 
collapse whereby particle states collapse into spacetime states. The extra feature here 
is the thermalization of the particle states produced through Hawking radiation, which 
eventually leads to the Information Loss conundrum. 

The issue we address here is that of thermal stability of generic radiant black holes 
in the heat bath made up of their own radiation. The asymptotically fiat Schwarzschild 
spacetime is well-known ^ to have a thermal instability: the Hawking temperature for a 
Schwarzschild black hole of mass M is given by T ~ 1/M which implies that the specific 
heat C = OM/dT < ! The instability is attributed, within a standard canonical ensemble 
approach, to the superexponential growth of the density of states p{M) ~ exp M ^ which 
results in the canonical partition function diverging for large M. 

The problems with an approach based on an equilibrium canonical ensemble do not ex- 
ist, at least for isolated spherically symmetric black holes, formulated as isolated horizons 
[8\ of fixed horizon area; these can be consistently described in terms of an equilibrium 
microcanonical ensemble with fixed A (and hence disallowing thermal fluctuations of the 



energy M). For A » Ipi^ncky i^ ^^^ been shown using Loop Quantum Gravity P, that all 
spherically symmetric four dimensional isolated horizons possess a microcanonical entropy 
obeying the Bekenstein-Hawking Area Law (BHAL) [Z],|3]. Further, the microcanonical 
entropy has corrections to the BHAL due to quantum spacetime fluctuations at fixed 
horizon area. These arise, in the limit of large A, as an infinite series in inverse powers 
of horizon area beginning with a term logarithmic in the area ^U] , with completely finite 
coefficients, 

3 

Smc = Sbh - - log Sbh + const. + C(S'^^) • (4) 

where Sbh = ^/4/pi„„cfc 

On the other hand, asymptotically anti-de Sitter (adS) black holes with spherical sym- 
metry are known jB] to be describable in terms of an equilibrium canonical ensemble, so 
long as the cosmological constant is large in magnitude. For this range of black hole pa- 
rameters, to leading order in A the canonical entropy obeys the BHAL. As the magnitude 
of the cosmological constant is reduced, one approaches the so-called Hawking-Page phase 
transition to a 'phase' which exhibits the same thermal instability as mentioned above. 

In this article, we focus on the following 

• Is an understanding of the foregoing features of black hole entropy and thermal 
stability on some sort of a 'unified' basis possible ? We shall argue, following [12]- 
[T7j that it is indeed so, at least in the case of non-rotating black holes. 

• In addition to corrections (to the area law) due to fixed area quantum spacetime 
fiuctuations computed using a microcanonical approach, can one compute correc- 
tions due to thermal fiuctuations of horizon area within the canonical ensemble ? 
Once again, the answer is in the affirmative. The result found in fn]-|IZj; at least 
for the leading log area corrections, turns out to be universal in the sense that, just 
like the BHAL, it holds for all black holes independent of their parameters. 

2 Canonical partition function : holography ? 

Following jT^I, we start with the canonical partition function in the quantum case 

Zc{p) = Tr exp-(3H . (5) 

Recall that in classical general relativity in the Hamiltonian formulation, the bulk Hamil- 
tonian is a first class constraint, so that the entire Hamiltonian consists of the boundary 
contribution Hs on the constraint surface. In the quantum domain, the Hamiltonian 
operator can be written as 

H = Hv + Hs, (6) 



with the subscripts V and S signifying bulk and boundary terms respectively. The Hamil- 
tonian constraint is then implemented by requiring 

Hv \i^)v = (7) 

for every physical state \ip)v in the bulk. Choose as basis for the Hamiltonian in (jHl) the 
states \ip)v ® \x)s- This implies that the partition function may be factorized as 

Zc = Tr exp —/SH 

= dim Tiv Trs exp —(3Hs (8) 

indep. of f3. boundary 

Thus, the relevance of the bulk physics seems rather limited due to the constraint ^. 
The partition function further reduces to 

Zci/3) = dim Hv ZsiP) , (9) 

where TCy is the space of bulk states \i/j) and Zs is the 'boundary' partition function given 
by 

ZsiP) = Trs exp-pHs. (10) 

Since we are considering situations where, in addition to the boundary at asymptopia, 
there is also an inner boundary at the black hole horizon, quantum fluctuations of this 
boundary lead to black hole thermodynamics. The factorization in eq.® manifests in 
the canonical entropy as the appearance of an additive constant proportional to dim, Tiy- 
Since thermodynamic entropy is defined only upto an additive constant, we may argue 
that the bulk states do not play any role in black hole thermodynamics. This may be 
thought of as the origin of a weaker version of the holographic hypothesis ^Hl • 

For our purpose, it is more convenient to rewrite (jiup as 

Zc{(^) = J29iEsiAin))) exp-(3Es{A{n)), (11) 

ne2^^ • ' 

degeneracy 

where, we have made the assumptions that (a) the energy is a function of the area of the 
horizon A and (b) this area is quantized. The first assumption (a) basically originates in 
the idea in the last paragraph of that black hole thermodynamics ensues solely from the 
boundary states whose energy ought to be a function of some property of the boundary like 
area. The second assumption (b) is actually explicitly provable in theories like NCQGR 
as we now briefly digress to explain. 




Figure 2: Internal boundary (horizon) pierced by spinet links 

3 Spin network basis in NCQGR 

The basic canonical degrees of freedom in NCQGR are holonomies of a distributional 
SU{2) connection and fluxes of the densitized triad conjugate to this connection. The 
Gauss law (local SU{2) invariance) and momentum (spatial diffeomorphism) constraints 
are realized as self-adjoint operators constructed out of these variables. States annihilated 
by these constraint operators span the kinematical Hilbert space. Particularly convenient 
bases for this kinematical Hilbert space are the spin network bases. In any of these bases, 
a ('spinet') state is described in terms of links /i, . . . , /„ carrying spins {SU{2) irreducible 
representations) ji,...jn and vertices carrying invariant SU{2) tensors ('intertwiners'). 
A particularly important property of such bases is that geometrical observables like area 
operator are diagonal in this basis with discrete spectrum. An internal boundary of a 
spacetime like a horizon appears in this kinematical description as a punctured S^ with 
each puncture having a deficit angle 6 = 6{ji),i = 1, . . . ,p, as shown in Fig.l. 

For macroscopically large boundary areas A » Ipianck^ the area spectrum is domi- 
nated by ji = 1/2, Wi = 1, . . . ,p, p » 1. This is the situation when the deficit angles 
at each puncture takes its smallest nontrivial value, so that a classical horizon emerges. 
That implies that 

Mp) ~ P lllanck , P » 1 • (12) 

This completes our digression on NCQGR. 

4 Fluctuation effects on canonical entropy 

We now move on to discuss the effect of inclusion of Gaussian thermal fluctuations of the 
horizon area. The canonical entropy is expected to receive additional corrections due to 
such fluctuations over and above those due to quantum spacetime fluctuations already 



included in the microcanoncal entropy. Going back to eq. (jllj) . we can now rewrite the 
partition function as an integral, using the Poisson resummation formula 



E /( 
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dx exp{—2TTimx) f{x) 



(13) 



For macroscopically large horizon areas A{p), a; >> 1, so that the summation on the rhs 
of (fTSj) is dominated by the contribution of the m = term. In this approximation, we 
have 
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dx g{E{A{x))) exp-pE{A{x)) 



dE exp[S'A/c(^) 



log|-| - PE] 



(14) 



where Sue = ^ogg{E) is the microcanonical entropy. 

Now, in equilibrium statistical mechanics, there is an inherent ambiguity in the def- 
inition of the microcanonical entropy, since it may also be defined as Smc = ^ogp{E) 
where p{E) is the density of states. The relation between these two definitions involves 
the 'Jacobian' factor \dE/dx\~^ 
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MC 



Smc - log 



dE 
dx 



(15) 



Clearly, this ambiguity is irrelevant if all one is interested in is the leading order BHAL. 
However, if one is interested in logarithmic corrections to BHAL as we are, this difference 
is crucial and must be taken into account. 

We next proceed to evaluate the partition function in eq. (|Tl|l using the saddle point 
approximation around the point E = M where M is to be identified with the (classical) 
mass of the boundary (horizon). Integrating over the Gaussian fluctuations around the 
saddle point, and dropping higher order terms, we get. 
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exp {Smc{M) -(3M- log \-^\e=m]} 



TT 



nl/2 



Using Si 
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log Zc + PM, we obtain for the canonical entropy Sc 
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Smc{M) - -log(A) 



(16) 
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where. 
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[A{x)f 



SmcK"^) 
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M'{A) 
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Thus, the canonical entropy is expressed in terms of the microcanonical entropy for an 
average large horizon area, and the mass which is also a function of the area. Clearly, 
stable equilibrium ensues so long as A > 0. 

Additional support for this condition can be gleaned by considering the thermal ca- 
pacity of the system, using the standard relation 



C{A)^ 



dM _ M'{A) 

It ~ T'{A) ' 



(19) 



with T being derived from the microcanonical entropy Smc{A), and hence a function of 
A. One obtains for the heat capacity the relation 



CiA) 



M'{A) 
T{A)A'{x) 



n2 



A- 



(20) 



so that C > if only if A > 0. Since the positivity of the heat capacity is certainly 
a necessary condition for stable thermal equilibrium, it is gratifying that an identical 
criterion emerges for A as found from the canonical entropy (fTTjl . 

Using now eq, (fTHj) for the expression for A, the criterion for thermal stability of 
non-rotating macroscopic black holes is then easily seen to be 



M{A) > Smc{A) 



(21) 



as already mentioned in the summary. We have been using units in which G = h = c = 
/cfi = 1. If we revert back to units where these constants are not set to unity, the lower 
bound eq. (PT|) can be re-expressed as 



M{A) > 



he 



1/2 



Smc{A) 



(22) 



We remind the reader that in contrast to semiclassical approaches based on specific prop- 
erties of classical metrics, our approach incorporates crucially the microcanonical entropy 
generated by quantum spacetime fluctuations that leave the horizon area constant. Apart 
from the plausible assumption of the black hole mass being dependent only on the horizon 
area, no other assumption has been made to arrive at the result. Even so, it subsumes 
most results based on the semiclassical approach. 

As a byproduct of the above analysis, the canonical entropy for stable black holes can 
be expressed in terms of the Bekenstein-Hawking entropy Sbh as 



S, 
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S',,c{A) M"{A) 
Sl,c{A)M'{A) 



(23) 
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For any smooth M{A), one can truncate its power series expansion in A at some large 
order and show that the quantity in square brackets in eq. (PH|) does not contribute to 
the log(area) term, so that 

Sc = Sbh -^ (^-1) ^^^Sbh + const. + OiS^h) , (24) 

where, ^ = 3 in eq. ^. Note that this is the result for an isolated horizon described by 
an SU{2) Chern Simons theory. For a U{1) Chern Simons theory, i^ = 1 JI], JHl- The 
interplay between constant area quantum spacetime fluctuations and thermal fluctuations 
is obvious in the coefficient of the log{area) term where the contribution due to each ap- 
pears with a specific sign. It is not surprising that the thermal fluctuation contribution 
increases the canonical entropy. The cancellation occurring for horizons on which a resid- 
ual U{1) subgroup of SU{2) survives, because of additional gauge fixing by the boundary 
conditions describing an isolated horizon [5^, may indicate a possible non-renormalization 
theorem, although no special symmetry like supersymmetry has been employed anywhere 
above. It is thus generic for all non-rotating black holes, including those with electric or 
dilatonic charge. 

While so far we have restricted our attention to thermal fluctuations of area due 
to energy fluctuations alone, the stability criterion (J2H) can be shown to hold when in 
addition thermal fluctuations of electric charge are incorporated within a grand canonical 
ensemble ^Hl- As in ^5, we assume that energy spectrum is a function of the discrete 
area spectrum (well-known in LQG J2]) and a discrete charge spectrum. The charge 
spectrum is of course equally spaced in general; for large macroscopic black holes the area 
spectrum is equally spaced as well. 

In a basis in which both the area and charge operators are simultaneously diagonal, 
the grand canonical partition function can be expressed as 

Zg = Y. 9{'m,n) exp -/5 [ E(A^, Q„) - $ Q„] , (25) 

•m,n 

where, g [171,11) is the degeneracy corresponding to the area eigenvalue Am and charge 
eigenvalue Qn- Using a generalization of the Poisson resummation formula 

Y. fi^^^) = Y. I dx dy exp {-i{kx + ly)} f{x,y) (26) 

m,n k,l 

and assuming that the partition sum is dominated by the large eigenvalues A^ ,Qn, it 
can be expressed as a double integral 

Zg = J dxdy exp - {3 {E{A{x),Q{y)) - ^ Q{y)} g{A{x),Q{y)) . (27) 

Note that the transition from the discrete sum to the integral for Zg requires only that 
the dominant eigenvalues are large compared to the fundamental units of discreteness 
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which for the area is the Planck area and for the charge is the electronic charge. These 
conditions are of course fulfilled for all astrophysical black holes. 
Changing variables in eq. ()27p form x,y to E,Q 



>G 



j dE dQ J{E,Q) g{E,Q) exp {-/5(E - <l>g)} 
jdEdQp{E,Q) exp {-/3(E - $g)} , 



(2^ 



where, the Jacobian J = \ E^^ \ \ Q,y \ , and p = Jl{E, Q) g{E, Q) is the density 
of states. Employing the saddle point approximation and using eq. (|H^ one obtains 



Sg{M,Qo) = SMciM,Qo) 



- log A + const. 



(29) 



where, using Smc = log P = Smc + log J we have defined 

= detfi(E,)2(gj2 
= deiVt {A^E? {A^^f {Q^yf\M,Q, 



A = detfijT^ = detfi(E^)2(g ^2 



y) \M,Qo 



where, the Hessian matrix 
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(30) 



(31) 



Since the micro canonical entropy is known to be only a function of the horizon area 
even for charged non-rotating black holes 1111101 , one can express det fl as 



detf] 



(SmcaY (Aee) (A^qq) - {A^eqY 



\M,Qo 



+ s 



MC,AA JMC,A 



Sa 



A,EE — 2 A^E A,Q ^,EQ 



(A^eY A,qq + (A 

The necessary and sufficient conditions for thermal stability are 

Trn -- 
det(] : 



\M,Q, 



o(32) 



Spic,ee\m,Qo + Smc,qq\m,Qo ^ ^ 



Sa 



Sa 



- s^ 



M,Qo 



niCEE OMCQQ - ^MC,EQ 

which necessarily imply 

Smc,ee\m,Qo < ) and Smc,qq\m,Qo < . 



> 



(33) 



(34) 



Note that while these conditions together imply the first of the necessary and sufficient 
conditions (jHHI) for stability of fl, they are not sufficient to guarantee the second one. 
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Using the microcanonical relations for temperature and potential, we may expess the 
necessary conditions for stability in terms of the heat capacity Cq = {dE/dT)Q and the 
capacitance C = {dQ/d^)E in the following way 

Cq > and C $ < T (^] . (35) 

The more stringent necessary and sufficient conditions can also be similarly expressed in 
terms of Cq and C. 

The grand partition function, evaluated in the saddle point approximation, can now 
be substituted in the standard thermodynamic relation in the presence of a chemical 
(electrostatic) potential 

Sg = PM-(3Q<!>+ logZc, (36) 

so as to yield the grand canonical entropy 

Sg = Smc — - log det Q + const. (37) 

We now make use of eq. (JH) to observe that both Smc,a and Smc,aa are positive 
definite for macroscopically large areas. The necessary and sufficient conditions for ther- 
mal stability, under both energy and charge fluctuations, can now be transcribed into the 
three simple inequalities 

^MM Smc,AA p, 

Am Smc,a 

Aqq Smc,aa q 



Art O 



MC,A 



2 



^MM AqQ ( Smc,AA \ /oo\ 

^M Aq \ Smc,a J 

where, we have used the notation A^ee\m,Qo — Amm etc. We have assumed of course 
that the horizon area A = A{M, Q). Since we are identifying the black hole mass M with 
the mass deflned for the isolated non-radiating non-accreting horizon, we can recall our 
earlier assumption that M = M{A), so that A = A{M{A),Q). In other words, the mass 
associated with the isolated horizon ought to relate to the horizon area just as the bulk 
Hamiltonian relates to the volume operator (as shown in ref. fH])- The quantum mass 
spectrum thus should be related to the area spectrum. On the other hand, electric charge 
has no such geometric origin, and is independent of the area, as far as one can make out. 
This implies that 

Am M'{A) = 1 
Amm M"{A) + Am M"{A) = . (39) 
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Substituting eq. (j39|) into the inequality (|38p. the stability criterion eq. (|2H) emerges 

once again. It is not difficult to see that this, when substituted back into (j38|) . merely 
reproduces (jHHj) . thereby establishing consistency between the inequalities. 

Is it possible to derive as a bonus, as in the case of non-fluctuating charge, a general 
formula for the thermal fluctuation correction to the canonical entropy, logarithmic in the 
horizon area ? We believe it is indeed the case, but do not include that discussion here. 

Can this criterion continue to hold if the black hole is characterized by a number of 
[/(I) charges Qi, . . . , Qn-i, all independently quantized ? This situation typically ensues 
in black holes arising in the low energy supergravity limit of string theories. The matrix Q 
is then an n x n matrix, and the analogues of eq.s (jnHj) -()38 p now become more complicated, 
involving sums of products of derivatives of Smc- Despite this, so long as one stays away 
from extremality, it is not inconceiveable that the mass alone decides on the stability 
issue. 

5 Concluding remarks 

The laws of black hole mechanics derived from the Raychaudhuri equation have inevitably 
led to worldwide attempts to seek quantum formulations of general relativity. A reason- 
ably satisfactory understanding of black hole entropy has been achieved : the entropy 
of radiant black holes is completely described in terms of the microcanonical entropy of 
isolated horizons which, in turn, is more or less entirely understood within the scheme of 
Loop Quantum Gravity, at least for the cases without rotation. Inclusion of rotation into 
the LQG approach to microcanonical entropy remains now foremost on the agenda, since, 
like the mass of the isolated horizon, and unlike f/(l) charges, the angular momentum is 
also likely to be determined by the horizon area in the isolated case. 

Acknowledgement 1 thank A. Chatterjee for several illuminating discussions. 
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